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The computation of mode sums of the types encountered in basic quantum field theoretic 
applications is addressed with an emphasis on their expansions into functions of distance 
that can be interpreted as potentials. We show how to regularize and calculate the Casimir 
energy for the continuum Nambu-Goto string with massive ends as well as for the discrete 
Isgur-Paton non-relativistic string with massive ends. As an additional example, we examine 
the effect on the interquark potential of a constant Kalb-Ramond field strength interacting 
with a QCD string. 



I. INTRODUCTION 



The determination of the Casimir, or zero-point, energy for a quantum system through the evaluation 
of the zero mode energy sum is a problem of quite practical importance. Casimir originally considered the 
effect of conducting plates on the QED vacuum energy due to the changed boundary conditions Q . A review 
of the Casimir effect can be found in the Physics Report by Plunien, Muller, and Greiner Q. 

The zero-point energy sums are almost always divergent and need regularization in order to be interpreted 
sensibly. The resulting Casimir energy may be interpreted as a background energy or potential. This is 
made most evident if the sum can be expressed in terms of simple functions of distance, which might be 
put into one-to-one correspondence with simple potential terms such as Coulombic or centripetal potential 
components. The purpose of this paper is to illustrate how this can be done for several classes of mode sums 
encountered recently in the literature. 

We begin in section by examining the computation of the Casimir energy for a static string with massive 
ends that is allowed to vibrate transversely. This problem was solved in closed integral form by Lambiase 
and Nesterenko ||. We review the derivation of the integral form and then find an approximate series 
representation for large quark end masses. 

In section II] , we find useful series approximation for a sum of Bessel functions that appears throughout 
this paper, and elsewhere in physics literature. 

In section IV we examine the effect of a particular background Kalb-Ramond field upon the Casimir 
energy. The Kalb-Ramond field we examine is one giving rise to a constant three-form field strength. The 
Kalb-Ramond field strength may be viewed as the dual of the abelian monopole current found in lattice QCD 
that likely leads to the formation of the chromoelectric flux tubes. The effect on the Casimir energy of string 
self- interaction mediated by the Kalb-Ramond field and the dilaton is to renormalize the string tension 
which also happens in a background Kalb-Ramond field. We find additional non-linear corrections to the 
static potential in a background Kalb-Ramond field. 

Finally, we examine the Casimir energies in the discretized string model of Isgur and Paton [T^ ] . We find 
the interesting result that although the Casimir energies of the continuum and discrete strings agree for 
fixed-end strings, the finite quark mass corrections are of a completely different form. 



II. CONTINUUM QCD STRING WITH MASSIVE ENDS 

The eigenfrequencies for a vibrating Nambu-Goto string of length r and tension a, terminated by two 
massive quarks of mass M that are allowed to move only transversely has been given by Lambiase and 
Nesterenko H as the roots of the transcendental equation 

f{uj) = (M 2 uj 2 - a 2 )sh\ujr - 2Maucosur = . (1) 

The Casimir energy, which is half the sum of these eigenfrequencies, can be calculated with aid of the 
Watson-Sommerfcld transform 

-^Casimir = ^ = ~7~~T / U duj In f(ij) , (2) 

2 z ^ 4tti J c duj 

where the contour C encloses the roots of f(u>), all of which lie on the positive real axis for physically 
meaningful frequencies. The energy computed in Eq. (|^) is the energy per transverse degree of freedom. In 
D dimensions, there would be an extra factor of D — 2. The sum in Eq. (Q) is formally infinite and must be 
regularized and renormalized. We deform the contour to include the entire left half-plane and keep only the 
regular portion 
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We renormalize by subtracting off the part of this integral for which r — > oo, to arrive at 
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and integrate by parts to arrive at 
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Defining the quantity i~(s), we have the exact result 
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The evaluation of the Casimir energy for small s can be performed as follows. For very small s (large masses 
M) the approximation, 
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is extremely accurate. The logarithm can be expanded 
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and used with the integral representation of the modified Bessel function 
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as an excellent approximation for s small. Using the results of section III below, we find 

7(s) = - — + V2tt 2 s + (47 - 2)s - 4s sintT 1 f -^=) + • • • , 
12 \V8s/ 

which agrees extremely well with numerical evaluations of the integral for small s, as can be seen in Fig. [l]. 
A comparison of the series ([l5]) and the sum ([TJ|) is shown in Fig. If more accuracy is needed for larger 
values of s, higher terms in s in the exponent of Eq. ( pTj| ) can be retained and those contributions to the 
exponent expanded in a power series. 

If the quark mass M is large, we can use the approximation Eq. ([l5|), to calculate the lowest order 
corrections in jj to the Casimir energy. From the definition Eq. (^), we find the Casimir energy per 
transverse degree of freedom 
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III. EVALUATION OF A SUM OF BESSEL FUNCTIONS 



In this section we evaluate the sum 

G(x) = V -K^nx) , (17) 
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for small a;. This sum appears in the evaluation of functional determinants Q as well as Casimir energies 
Hf||7j . Related sums of Bessel functions appear in generic one-loop finite temperature effective potentials 
for GUTs H , and also appear in the mathematics literature 

Evaluation of this sum using zeta functions quickly becomes problematic because no matter how small x 
is, for sufficiently large n values, nx — > oo and small x expansions of K\(x) cannot be used. However, we 
can use another integral representation of the Bessel function, 
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together with the sum 
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where [t] is the greatest integer less than t and B2(t) is the second Bernoulli polynomial, 

B 2 {t) = \ - t + t 2 . (20) 
o 

We obtain 
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We may split this expression into the first term and a remainder 

G(x) = 2., T ^ ,* + 2,x r ^-^\dt . (22) 
Jo {x 2 + (2vt) 2 )i ii (a; 2 + (27rt) 2 )i 



The remainder may be written as a sum, 
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which may then be expanded for small i in a binomial series; 

2wx y f 1 _!____> dt-2^xyH)-^—y f 1 ____! dt (24) 

Finally, we rewrite this expression in term of the generalized zeta function 

G(x) = 2nx 2 A±^ dt + 2iT 3 xy 2 / B 2 (t)((3 + 2k, t) dt . (25) 
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All of these integrals are strongly convergent, which is obvious from simple power counting for all but the 
k = term, which we will explicitly evaluate. The integral in the first term, which is dominant for x w 0, is 
elementary 
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For the fc = term, we replace the sum of integrals by its original form to obtain 
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We evaluate each of the integrals separately. The first integral is trivial and gives Using the notation 
F-zim) = m~ 2 , we find that the second integral is 
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We can evaluate the last integral similarly, this time using the notation F\ (m) = m 

— ' m dt 



/ ^ dt = lim lniV-2y / rr+> / 

A t 3 N^oo\ Z^J Q {t + m) 2 Z^J Q 



Af^ooji t 3 N^oo \ ^ Jn (t + m) 2 ! Jn (t + m) 3 

\ m—1 u v 7 m—1 u v 7 / 

CAT— 1 _^ JV-1 \ 

lniV + 2 m(Fi{m + 1) - Fi(m)) - - ^ m 2 (F 2 (m + 1) - F 2 (m)) j 
m — 1 m—1 / 



. — (30) 

m—1 / 

Rearranging the series according to Eq. (|28|) and 
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We put all the terms together from Eqs. (|27|), (|29|), and fl33| ) to find 
which implies that 
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IV. QCD STRING INTERACTING WITH A BACKGROUND MONOPOLE DENSITY 



The physics of charged particles moving in a plane under the influence of a constant magnetic field is a 
very beautiful and enlightening subject containing such surprises as the quantum Hall effect. The string 
analog of a charge moving in a constant magnetic field is the motion of a string in a constant Kalb-Ramond 
g§ field strength 

H 0ij (x) = , 

Hijk(x) = heij k . (36) 
This system is a model of a QCD string interacting with a constant monopole density and, like the Landau 



problem, is also exactly solvable |p.l| . (The monopole current is proportional to the dual of i? Q /3 7 : J, 
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e A ' a ^ 7 .H Q , | g 7 .) The Kalb-Ramond potential B^ y corresponding to the field strength (J3g) is 

B 0l {x) = , 

B ij(x) = ^jkx k . (37) 

In vortex mechanics, the Kalb-Ramond field describes the Magnus forces on a vortex line moving in a fluid. 
The first-quantized action for a bosonic string moving in a background Kalb-Ramond field is 



| J ' dadr (^hh ab d a X»d b X v Vllv + B^(X)e ab d a X»d b X v ) . (3N) 



If we substitute the field Bij above into the string action, and we choose orthonormal coordinates on the 
worldsheet, we find the action becomes 

S=-^Jdadr (v ab d a X»d b X v + K l3k e ab X l d a X^d b X k ^j , (39) 

and implies the equations of motion 

X° - X°" = , (40) 

X - X" = -KX x X' . (41) 

Here we have used reparametrization invariance of the action to set coordinate time equal to laboratory 
time t — x° and to set the spatial parameter of the string equal to the distance in real space along the 
string, a — z. We wish to find the Casimir energy of a long straight string acting under the influence of this 
external Kalb-Ramond field. We assume the string lies along the z direction and has length r. Thus we find 
that X' = z, and the linearized equations of motion for the small vibrations on the string become 

X - X" = -hY , 

Y-Y" = hX . (42) 

We define the right and left circularly polarized modes X ± = X ± iY, and find that in terms of them the 
small amplitude vibration equations above reduce to 

X± _ x ±/i = Tih x± _ (43) 

Since our string has fixed length, we know that the wavelength with either fixed Dirichlet or Neumann 
boundary conditions (infinitely massive or completely massless quarks at the end) lead to 

X ± ± ihX ± = X ±n = - [— J X ± , n = 1,2,3,... . (44) 
or, if we assume that X = e lu>t Xo, 

c4±/^ ± =(^) 2 = A£. (45) 
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The frequencies of small vibrations are then 
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Naively, the Casimir energy as a function of the Kalb-Ramond field strength h is half the sum of all the 
frequencies 
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Mode sums of this form, 
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appear in a large variety of problems (see the references following Eq. (5.33) in H). We now show the simple 
steps with which this sum can be expressed in a form very similar to Eq. (14). 
First, we subtract out the divergence from the sum in order to regularize it 
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The sum may be done by Poisson resummation 
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Concentrating on the last term, we may use 
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where K u is a modified Bessel function. Letting s = tx 2 , we rewrite the last term in the sum as 
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Finally, we have 
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in which the final integral may be expanded and evaluated in terms of Gamma functions as 
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Subtracting the infinite pieces above, we find the final regularized result 
a^(0 -- 
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This and related results are found by a slightly different method in Ref. . 
The Casimir energy corresponding to the zero-point frequencies ((46|), 
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can be evaluated using the results of section We find 

+ 0(h 3 r 2 ) . 

For h — > this reduces to the fixed-end straight string result for two transverse degrees of freedom. 
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V. DISCRETE ISGUR-PATON STRING 

In this section we consider the interquark potential of a static quark-antiquark pair bound by a flux tube, 
considered by Isgur and Paton |H| and Merlin and Paton jl3j to be a discrete string consisting of N masses m 
a distance d apart. This string, when terminated by particles of mass M, has equations of motion JT^ , p^| , ^4[ 
with the secular equation 

M sin(jVg) 

m sin((JV + l)0) ' [ 1 

uj 2 

A = 2(1- cos 0) , A = — , (62) 

P 

for the transverse vibrational frequencies uj. Here (3 is defined in terms of the masses m, the bare tension ao 
and interquark distance r, as 8 = g0 ( jy+1 ) . xhe infinite mass M limit results in 
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Assuming finite M corrections of the form 



p=l,2,...,iV. (63) 



(FTT) ' (64) 



7 



we discover that 
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The contribution of the first and third terms to the mode sum is simple, following from 
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The second term can be analyzed by use of the expansion 
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in which _B/j are the Bernoulli numbers. 
We begin with 
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in which the Bernoulli polynomials B n (x) satisfy the formal recursion 

B n {x) = (B + x) n with B n -> B n . 

We find that 
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The second term in Eq. ( |65| ) thus contributes five terms to the mode sum 
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of leading orders in N. The first is simple 
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for which A2 = A3 = A4 = ||, and so forth. This is a standard expansion of the Digamma function [p| 
The next term is 
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and subsequent terms can all be expressed as valuations of the derivatives of csc(x) — — at the point x = 
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and finally 
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All of the parts can be assembled to give 
J(N) 
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This expression agrees very well with the actual sum, plotted as a function of N. This is illustrated in Fig. ^. 

For large quark masses M, we can again calculate the lowest order i corrections to the Casimir energy, 
this time for the discrete string before the limit of link spacing going to zero (d — > 0) is taken. To make 
contact with Isgur and Paton fl2|| , we take f) = and then use AT=(^) — lto arrive at the Casimir energy 
contribution per transverse degree of freedom 
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We note that although the universal Luscher term — is present in both (|l6| ) and (|S2"|), the A corrections 
are completely different. In the former case, the correction is 0{-^=), while in the latter it is 0(-^). Of 
course, there is no reason to expect that the form of the Casimir energies will be the same for the Nambu- 
Goto string and the discrete string of Isgur and Paton. The former is relativistic while the later is not. In 
addition, it is not likely that the Casimir energy of a discrete string will go the the continuum result as the 
number of mass points is taken to infinity. Sums and limits generally depend upon the order in which they 
are done, especially when the sums are divergent and must be regularized. 
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VI. CONCLUSIONS 



We have shown how to regularize and perform the zero-point energy mode sums for the discrete non- 
relativistic string studied by Isgur and Paton and the relativistic continuum Nambu-Goto string. Computa- 
tion of these divergent sums is necessary in order to obtain the Casimir energy of the system to be interpreted 
as part of the static potential for quarks bound by a QCD string. We have shown that the lowest order 
corrections to the Casimir energies in the reciprocal of the mass of the terminating quark is different for the 
two cases. We have also shown how the addition of a background Kalb-Ramond field with a corresponding 
constant field strength affects the static potential. We hope to extend these results to the case of a string 
interacting with itself through a Kalb-Ramond field. 
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FIG. 3. Numerical evaluations of the sum J(N) . The solid curve is the function Eq. (glj), while the diamonds are 
for the sum Eq. (|^). 
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